Atom in a coherently controlled squeezed vacuum 
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A broadband squeezed vacuum photon field is characterized by a complex squeezing function. 
We show that by controlling the wavelength dependence of its phase it is possible to change the 
dynamics of the atomic polarization interacting with the squeezed vacuum. Such a phase modulation 
effectively produces a finite range temporal interaction kernel between the two quadratures of the 
atomic polarization yielding the change in the decay rates as well as the appearance of additional 
oscillation frequencies. We show that decay rates slower than the spontaneous decay rate can be 
achieved even for a squeezed bath in the classic regime. For linear and quadratic phase modulations 
the power spectrum of the scattered light exhibits narrowing of the central peak due to the modified 
decay rates. For strong phase modulations side lobes appear symmetrically around the central peak 
reflecting additional oscillation frequencies. 

PACS numbers: 78.67.De, 42.50.Dv, 42.55.Sa, 42.50.Lc 



The effect of the interaction of an atom with a squeezed 
light had been studied extensively. Gardiner |lj has con- 
sidered the behavior of a two-level atom damped by an 
infinite bandwidth squeezed vacuum. He showed that 
the two quadratures of the atomic polarization decay at 
different rates. Decay rates smaller than the decay rates 
of spontaneous emission can be achieved for non-classical 
squeezing. The presence of the two decay rates modify 
the fluorescence spectrum, Ref. Q. In Refs. [1 Q 
the interaction of a finite bandwidth squeezed vacuum 
with a two level atom was investigated. Review of these 
and related studies is found in Ref. Q. Recently these 
works were extended in Ref. Q to include interactions 
with semiconductor microstructures. 

In the field of quantum coherent control pulse shapers 
were used to attain prescribed phase modulation of a 
down converted light in order to control two photon ab- 
sorption. It was shown Q that pulses can be shaped in a 
way that will stretch them temporally affecting the tran- 
sition probability ^j- Experiments 0, El have also 
been performed on two photon absorption with coher- 
ent, narrow band down-converted light, demonstrating 
nonclassical features which appear at very low powers 
[3 Q E3 an d result from time and energy correla- 
tions (entanglement) between the down-converted pho- 
ton pairs. 

In this work we wish to investigate the effect of control- 
ling and modulating the relative phase of the modes of a 
squeezed reservoir interacting with an atom. This can be 
easily done by a pulse shaper arrangement as shown in 
Refs. We will demonstrate that by controlling in 

this manner the phases of the correlations in the squeezed 
reservoir it is possible to control the dynamics of atomic 
polarization and in particular to further reduce its decay 
rates. 

We will use the standard model to describe the in- 
teraction of a two-level atom with a broadband radiation 
field. The atom is assumed to be coupled to a one dimen- 



sional set of radiation modes. The Hamiltonian is given 
in the electric-dipolc and rotating-wave approximations 
by H = H + Hi as (we take % = 1): 

H = uj a a z + ^ u q b q b q 
i 

Hi = r<T+ + <T-ri , r = J2g q b q - (i) 

<? 

The pseudospins er + ,er_and a z describe the atom, uj a is 
the atomic resonance frequency, b q and b q describe radi- 
ation modes with wave vector q and frequency Lu q and g q 
are the mode-atom couplings. 

We assume that the radiation acts as a reservoir with 
correlations of a two mode broadband squeezed vac- 
uum dm 

{b\b q .) = N{u> q )8 q , q , 

(b\b\,) = M(u q )6 2qo - q<q , (2) 

(b q ) = (ftt) = o, 

where N(uj q ) describes the average occupations of pho- 
tonic modes while the magnitude |M(w g )| gives the 
squeezing strength of mode pairs centered around the 
frequency ujq = uj{qo)- The phase of M(uj q ) describes the 
"direction" of squeezing in the phase space of squeezed 
mode pairs. In the common method of generating 
squeezed vacuum radiation by non-linear down conver- 
sion the phase of M(uj q ) is constant over the frequency 
range. Letting the radiation pass through a pulse shaper 
makes it possible to change this phase into a prescribed 
function of uu q , cf., Refs. 

We shall for simplicity assume that N, \M\,g q are con- 
stants within a bandwidth ojq ± B/2 and that B is <C ojq 
but is much larger than any other frequency in the sys- 
tem. 
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Using the equations of motion 



and the following function 



-i[a+,H] 
-i[a z ,H] 
-i[b q ,H] : 



iu) a <J+ — 2iT^a z 
— iTa+ + iT^a- 



(3) 



-ILUqbg 



and transforming them to the rotating frame of the laser 
frequency, we integrate over time the equations for a z 
and b q and substitute the result back into the equa- 
tion for er+. We then average the resulting equation 
over the initial state and assume that the bath and 
the atom are (approximately) decorrelated, i.e. that 
the photon atom correlators factorize at all times e.g. 
(b q (t)a + (t')) w (b q (t))(a + (t')}, etc. The decorrelation as- 
sumption is valid in the regime of weak coupling between 
the system and the photon bath. Using it we obtain a 
closed equation for the atomic polarizations (<r±) 



df 



(M = (i(« B -W0)-|)<O-+> (4) 

- 2 f dt'{Y,\9 q \ 2 N(cUq)e^-^ ><*-*') (o+M) 

JO q 

where ujq = 2ujo—ujq, 7 = p(uj)\g q \ 2 is the vacuum atomic 
decay rate and p(u>) is the density of the radiation modes. 
We assume that p(oj) is flat over the bandwidth B. 

We now transform the sums over q into integrals. We 
define the following parameters 



1 N = P (u J )\gq\ 2 N(0Jq) 

1 M = p(uj)g* q g* Q \M(u q )\, 



(5) 



i-B/2 



k(t-t') = - f ' dwe^^e*"**-*'). (6) 

7T J-B/2 



where f(w) is the phase of M(u) = |M|e i/(w "" o) . Note 
that f(u>) — f(—u>). Without loss of generality we can 
take /(O) = using the freedom to absorb its non-zero 
value in the phase of M. . 

With this notation the equation for the atomic polar- 
izations becomes 



d_ 

dt 



(M = (iS-W+l)) (<t+)+^mJ dt'k(t-t'){a-(t')), 

(7) 

where we defined the detuning S = uj a — ujq . 



The phase modulation of the squeezing parameters 
leads to a finite range memory kernel in Eq. Q which 
couples atomic polarizations to their complex conjugates 
at earlier times. This non-Markovian dynamics is still lin- 
ear due to the atom-bath decoupling assumption. One 
can envisage situations in which the non-Markoian dy- 
namics may lead to the breakdown of this assumption 
but we will not consider such cases here. 

To analyze the controlled memory effects we apply the 
Laplace transform to Eq. Q and obtain 



(*+(*)} 



s + i6 + ~f(Af+ 1/2) 
jM*k*(s*) 



iS + y(J\f- 



1/2) 



(Mo)) 
(Mo)) 



(s + 1 {N + 1/2)) 2 + S 2 - \M\ 2 k(s)k*{s* 
I 



(8) 



where we assumed non-zero initial conditions at t = for 
(a±(t)} and denoted (cr±(s)), k(s) the Laplace transform 
of the polarizations and the kernel function respectively. 
The non-zero (c±(0)) is the easiest situation to analyze 
although it can only be realized with a specially designed 
initial pulse applied before letting the squeezed reservoir 
interact with the atom. We will later discuss the implica- 
tions for the fluorescence spectrum. In the following we 
will assume for simplicity the exact resonance case 8 = 0. 

The poles of the Laplace transform are the solutions 

of 



s/7 = -TV - 1/2 ± \M\yk(s) k*(s*) 



For the unmodulated phase f(to) — we obtain the usual 
Markov result k(t) — 8(t), k(s) = 1, in which case s± — 
7(yV+l/2±7|.M|), representing the splitting of the decay 
rate into the fast and the slow components under the 
influence of the squeezing, Ref. lj. For the modulated 
phase it is convenient to consider the graphic solution, 
cf. Fig^ The symmetry of /(w) implies that fc(0) = 1. 
Therefore for analytic k(s) to the first order in s, fc(s) « 
1 - s J °° tk{t)dt = 1 - ski. The right hand side of Eq.Q 
reduces to -TV - 1/2 ± \M\(1 — sRefci) and if Refei > 
then the slow component of the decay rate 



(9) 



s_ = -7(l-7|M|Refci)(TV + l/2-|M|) + ... 
can become even slower than in the Markov case. 



(10) 
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FIG. 1: Graphical solution of Eq. Q with N = M — 1 for 
the quadratic phase modulation, Eq. 11121 . The wavy lines 
show the right hand side of the equation for Tj — 5 (solid 
line) and T7 = 1 (dashed line) while the horizontal dotted 
lines are for T = - the standard Markov case. 



FIG. 2: The fluorescence spectrum for the quadratic phase 
modulation, Eq. 11 U and Af = M. = 1. Solid, dashed and 
dotted lines are for T7 = 5, T7 = 1 and T — (the Markov 
case) respectively. 



We wish to remark that in addition to the real poles, 
the complex poles of Eq. JSJ in the Res < half plane 
also play a role in determining the dynamics of the atomic 
polarizations. Their imaginary parts represent additional 
oscillation frequencies. If their real part is smaller than 
the real poles they may dominate the long time decay. 
Finally possible cuts may also make a contribution. Wc 
will discuss such contributions in the examples below. 

We consider two concrete examples of the phase mod- 
ulation. The simplest case to calculate is the quadratic 
phase 



2, .2 



f(u) = T z u 



(11) 



where T is a real positive parameter. The resulting 
Laplace transformed memory kernel in the B — ► 00 limit 
is 



k(s) 



erfc(V-iTs) 



(12) 



The solution of 10 with this k(s) is shown graphically by 
plotting the two sides of this equation in Fig^ One sees 
the increase of the fast and the decrease of the slow decay 
rates caused by the phase modulation. In this example 

7Re/ci in ^ is y^f Tj. 

For k(s) given by Eq. (|El|) the structure of Eq. © 
is much richer than that implied by the real poles. One 
finds in addition numerous complex poles, the positions 
of which depend on the parameters M, Ai, T and 7. 
Their influence can be seen in Fig. [2] where we plot 
(<r_(s = ioj)) I (<T-(s = iu>o))- Together with the signifi- 
cant narrowing of the central part of the peak relative to 
the Markov case T = one observes the appearance of 
the side lobes which reflect a complicated pole structure 
in the Res < half plane. 

The quantity plotted in Fig. can be related to the 
fluorescence spectrum of light emitted by the atom into 
empty modes of the radiation field. This spectrum is 



given by 



S(u) = J-Re 



-(to 



-(to)), 



(13) 

where the subscript ss denotes steady state. This expres- 
sion involves the average of the two time product of po- 
larization operators (cr+(fo +T")c-(£o))- It is not difficult 
to show that for to = and io+T = t this average satisfies 
the same dynamical equation Q as a single time average 
(a+(t)) provided that one uses the atom - bath decorre- 
lation assumption. Irrespectively of the initial condition 
the system relaxes to a steady state regime. One can 
therefore replace the initial time, i.e. the low limit t = 
in the integral in Eq. (JJJ by t provided it is chosen af- 
ter the steady state is reached. The initial condition is 
then given by (a + (t )a_(t )) sa = 1/2+ (<T z (t )) ss . The 
Laplace transformed solution of this set is given by © 
with (d±(f-)(s) replacing a±(s) in the left hand side and 
(cr±(0)) replaced by 

<M0>>^M^)>,, = i-^, (14) 



(<7+(0)> (<7+(t )<7+(to)>* 



0. 



The fluorescence spectrum is given by S(u>) = 
^Jtef (iuS)} , the normalized form of which is the 

quantity plotted in Fig. [3 

We note that the discussion above essentially means 
that under the set of the adopted assumptions the quan- 
tum regression theorem can be applied despite the finite 
memory effects induced by the phase modulation. 

In Fig. G3we show how the phase modulation increases 
the sensitivity to squeezing. The graphs represent three 
possible reservoirs: squeezed phase modulated, squeezed 
unmodulated and white noise reservoir. One observes 
that the narrowing of the central peak is present even for 
the reduced values of M provided a strong phase modu- 
lation is applied. 

We will now briefly discuss another example, the linear 
phase modulation 



(15) 
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cence spectrum. We observe features similar to those 
found in the quadratic modulation case and in particular 
the narrowing of the spectrum around the central fre- 
quency and the development of the side lobes. From an 
experimental point of view the advantages of the linear 
modulation is the relative ease of achieving it in practice. 
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with Y. B. Lcvinson. 



FIG. 3: Sensitivity of the fluorescence spectrum to the degree 
of squeezing. T 7 = 2, M = 0.5, (solid line), T = 0, M = 0.5, 
(dashed) and M. = (dotted), white noise reservoir. For all 
graphs M = 1. 




FIG. 4: Linear phase modulation 115H for Af = 1, M. = 1 and 
T7 = 5 (solid line), Tj = 1 (dashed) T = 0, (doted), the 
normal Markov case. 



where T is a real positive parameter. In this case k(s) 
is multivalued with s = as a branch point. The 
pole structure is now more involved and should be dis- 
cussed together with the branch cuts in the complex 
Res < half plane. In particular the discussion fol- 
lowing Eq. 10 should be modified. We will not dis- 
cuss it here but rather present in Fig. 0] the quantity 
(er_(s = iu>)) / (&-(s — iloq)) f° r this phase modulation. 
As discussed above this quantity represents the fluores- 



[1 

[2 

[3. 

[5 
[6 

[-: 

[9 
[10 

[11 
[12. 

[13 
[14 

[is; 

[16 



* Electronic address: itay.rabinak@weizmann.ac.il 
C. W. Gardiner, Phys. Rev. Lett. 56, 1917 (1986). 
H. J. Carmichael, A. S. Lane, and D. F. Walls, Phys. 
Rev. Lett. 58, 2539 (1987). 

C. W. Gardiner and M. J. Collett, Phys. Rev. A 31, 3761 
(1985). 

C. W. Gardiner, A. S. Parkins, and M. J. Collett, J. Opt. 
Soc. Am. B 4, 1683 (1987). 

A. S. Parkins and C. W. Gardiner, Phys. Rev. A 40, 3796 

(1989) . 

D. J. Dalton, Z. Ficek and S. Swain, J. Mod. Opt. 46, 
379 (1999). 

E. Ginossar and S. Levit, Phy. Rev. B 72, 075333 (2005). 
M. M. Wefers and K. A. Nelson, Opt. Lett. 20, (1995). 
D. Meshulach and Y. Silberberg, Nature 396, (1998); 
Phys. Rev. A 60, 1287 (1999). 

N. Dudovich, B. Dayan, S. M. Gallagher Faeder, and Y. 
Silberberg, Phys. Rev. Lett. 86, 47 (2001). 
N. P. Georgiades et al., Phys. Rev. Lett. 75, 3426 (1995). 
N. P. Georgiades, E. S. Polzik, and H. J. Kimble, Phys. 
Rev. A 55, R1605 (1997). 

J. Gea-Banacloche, Phys. Rev. Lett. 62, 1603 (1989). 
J. Javanainen and P. L. Gould, Phys. Rev. A 41, 5088 

(1990) . 

H.-B. Fei et al, Phys. Rev. Lett. 78, 1679 (1997). 

B. Dayan, A. Pe'er, A. A. Friesem, and Y. Silberberg, 
Phy. Rev. Lett. 93, 023005 (2004). 



